luogu2805 [NOI2009]植物大战僵尸

发表于 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
选某个植物就要选保护它的所有植物
典型的最大权闭合子图
愉快的打完了
一测样例
225?
发现出现了一个环（良心样例不多见了）
怎么去掉环呢
tarjan?
发现环连出去的点都相当于无敌的
所以哪些点是可以被攻击的呢
搜索！
发现奇怪的性质。。
这是topsort!
愉快的拓扑一遍，然后重建边
再dinic跑一遍
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
#include<bits/stdc++.h>
#define mp make_pair
#define pb push_back
using namespace std;
typedef long long LL;
typedef pair<int,int> PII;
inline LL read()
{
	LL x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}
	return x*f;
}
const int N=608;
int n,m;
int v[N],d[N];
int x,y,z;
int s,t;
int ans;
inline int getid(int x,int y)
{
	return (x-1)*m+y;
}
int nume,head[N],cur[N],headx[N];
struct node{
	int to,nxt,f;
}e[N*N*4],ex[N*N];
inline void addedge(int x,int y,int z)
{
	e[++nume]=(node){y,head[x],z};head[x]=nume;
	e[++nume]=(node){x,head[y],0};head[y]=nume;
}
inline void insert(int x,int y)
{
	++d[y];
	ex[++nume]=(node){y,headx[x]};headx[x]=nume;
}
bool del[N];
int q[N],he,ta;
inline void topsort()
{
	for(int i=1;i<=n*m;++i){
		if(!d[i])	q[++ta]=i;
		else del[i]=1;
	}
	while(ta){
		int x=q[ta];del[x]=0;
		--ta;
		for(int i=headx[x];i;i=ex[i].nxt){
			--d[ex[i].to];
			if(!d[ex[i].to])	q[++ta]=ex[i].to;
		}
	}
}
inline void rebuild()
{
	nume=1;
	s=0;t=n*m+1;
	for(int i=1;i<=n*m;++i){
		if(!del[i]){
			if(v[i]>0){
				ans+=v[i];
				addedge(s,i,v[i]);
			}
			else{
				addedge(i,t,-v[i]);
			}
			for(int j=headx[i];j;j=ex[j].nxt){
				if(!del[ex[j].to]){
					addedge(ex[j].to,i,INT_MAX);
				}
			}
		}
	}
}
int dis[N];
inline bool bfs()
{
	memset(dis,0,sizeof(dis));
	he=1;ta=2;
	q[1]=s;
	dis[s]=1;
	while(he!=ta){
		int x=q[he];
		++he;
		for(int i=head[x];i;i=e[i].nxt){
			if(e[i].f&&!dis[e[i].to]){
				dis[e[i].to]=dis[x]+1;
				if(e[i].to==t)	return 1;
				q[ta]=e[i].to;
				++ta;
			}
		}
	}
	return 0;
}
inline int dfs(int x,int low)
{
	if(x==t||!low)	return low;
	int flow=0,tmp;
	for(int i=head[x];i;i=e[i].nxt){
		if(dis[e[i].to]==dis[x]+1&&(tmp=dfs(e[i].to,min(low,e[i].f)))){
			flow+=tmp;
			low-=tmp;
			e[i].f-=tmp;
			e[i^1].f+=tmp;
			if(!low)	return flow;
		}
	}
	if(low)	dis[x]=0;
	return flow;
}
inline int dinic()
{
	int maxflow=0;
	while(bfs()){
		for(int i=s;i<=t;++i){
			cur[i]=head[i];
		}
		maxflow+=dfs(s,INT_MAX);
	}
	return maxflow;
}
int main()
{
	n=read();m=read();
	for(int i=1;i<=n;++i){
		for(int j=1;j<=m;++j){
			v[getid(i,j)]=read();
			z=read();
			while(z--){
				x=read();y=read();
				++x;++y;
				insert(getid(i,j),getid(x,y));
			}
		}
	}
	for(int i=1;i<=n;++i){
		for(int j=2;j<=m;++j){
			insert(getid(i,j),getid(i,j-1));
		}
	}
	topsort();
	rebuild();
	ans-=dinic();
	printf("%d",ans);
	return 0;
}